Feshbach Resonance Induced Fano Interference in Photoassociation 
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We consider photoassociation from a state of two free atoms when the continuum state is close to 
a magnetic field induced Feshbach resonance and analyze Fano interference in photoassociation. We 
show that the minimum in photoassociation profiles characterized by the Fano asymmetry parameter 
q is independent of laser intensity, while the maximum explicitly depends on laser intensity. We 
further discuss the possibility of nonlinear Fano effect in photoassociation near a Feshbach resonance. 
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1. INTRODUCTION 

In recent times, quantum interferences have occupied a 
prominent place in physics and these occur rather ubiqui- 
tously. Many well-known examples of these include Fano 
interferences plfl E| , electromagnetically induced trans- 
parency (EIT) [J] jVacuum induced interferences in spon- 
taneous emission [5j] . The quantum interferences have re- 
sulted in large number of applications in coherent control 
of the optical properties, control of spontaneous emis- 
sion [6[ and slow light [Tj, [8( . Quantum interference has 
been experimentally demonstrated in coherent formation 
of molecules Q and Autler-Townes splitting [13, HH in 
two-photon PA. Theoretical formulation of PA within 
the framework of Fano's theory has been developed in 
Refs. [l2j and fl3l. Recent experimental [III [III, [ly] and 
theoretical jl7l. Il8l. Il9l . [2fj| studies on photoassociation 
(PA) near a magnetic field Feshbach resonance (MFR) 
[2lj have generated a lot of interest in Fano interfer- 
ence with ultracold atoms. In a remarkable experiment, 
Junker et al. [l4| have demonstrated asymmetric spectral 
line shape and saturation in PA due to a tunable MFR. 
Asymmetric line shape is a hallmark of Fano-effect and 
the experimental results of [l4[ can be attributed to the 
Fano interference. 

Here we demonstrate quantum interference in the con- 
text of photoassociation (PA) [22j under the condition 
when a Feshbach resonance is also involved in photoas- 
sociation. We show Fano like interference minimum in 
photoassociation spectrum. In analogy to the well know 
Fano q-parameter we can introduce a parameter which 
governs the existence of this minimum. Although the 
minimum is independent of laser intensity, the maximum 
is shown to depend explicitly on laser intensity. From our 
calculations we extract line shapes which are in broad 
agreement with the experimental results of junker et al. 
[p4| . Our formula for photoassociation is expressed in 
terms of parameters each of which has a clear physical 
meaning and is measurable. We derive probability of PA 
excitation for arbitrary intensities of the laser field and 
thus we also discuss nonlinear Fano effect. The current 




FIG. 1: A Schematic diagram showing the diatomic poten- 
tials, continuum and bound states, bound-bound and bound- 
continuum laser couplings. Note that the two double-arrow 
vertical lines refer to the same laser field - green line corre- 
sponds to continuum-bound and the magenta line refers to 
bound-bound couplings. The ground asymptotic channels 
(f a + fb) in the absence of a magnetic field corresponds to 
separated atoms in two different hyperfine numbers f a and 
fb- 



work has some features in common with the recent paper 
of Kuznetsova et al. [20( though these authors address a 
different problem which is the population transfer using 
two laser beams. Our emphasis is on quantum interfer- 
ences in PA using a single laser beam. 

The paper is organized in the following way. In sec- 
tion 2, we consider a simple model of three-channel time- 
independent scattering in the presence of an optical and 
a magnetic field. By using Green's functions, we present 
compact analytical solution of the model. We then dis- 
cuss selective results in Sec. 3. The paper is concluded in 
Sec.4. 



2. THE MODEL AND ITS SOLUTION 

To begin with, we model PA in the presence of a Fes- 
hbach resonance as a three-channel scattering problem. 
There are two ground-state asymptotic hyperfine chan- 
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nels of which one is closed and the other one is open. The 
third channel corresponds to the photoassociated excited 
molecular configuration. The two ground state channels 
are coupled via hyperfine interaction. At a Feshbach res- 
onance, the two atoms will form a quasibound state in 
the closed channel as schematically illustrated in Fig.l. 
As the strength of the applied magnetic field is varied, 
this quasibound state can move across the collision en- 
ergy. When a PA laser is applied to form an excited 
photoassociated molecule (PM), there arise two compet- 
ing pathways of dipole transitions as shown by different 
colors in Fig.l. One is the continuum-bound and the 
other one is bound-bound transition. We assume that 
the energy spacing of closed-channel quasibound states 
and rotational spacing of PM states are much larger than 
PA laser line width so that only one rotational level (J) 
of a particular vibrational state v of PM is coupled to a 
particular quasibound state by the PA laser. 

Let us write an energy eigenstate of the system of two 
atoms interacting simultaneously with a magnetic and a 
PA laser field in the form 



^ E ) = \ 92} +X \ 9i) + % \ e) 



(1) 



where E is an energy eigenvalue, | 51(2)) represents the 
internal electronic states of 1(2) or open(closed) chan- 
nel and I e) denotes the electronic state of the excited 
molecule. $/ and <& p are the diatomic bound states. The 
continuum state has the form \ — J dE'bE'ipE' where 
ipE' is an energy-normalized scattering state of collision 
energy E' and bE> is the density of unperturbed con- 
tinuum states. The state (fTJ) is assumed to be energy- 
normalized. The Hamiltonian of the system can be writ- 
ten as H = H kin + H dec + H hfs + H B + H L where 
Elkin denotes a term corresponding to the total kinetic 
energy of the two atoms and H e i ec is a term that de- 
pends on only electronic coordinates of the two atoms, 
Hhfs is the hyperfine interaction term. Here Hb repre- 
sents the magnetic interaction in the atomic states, and 
Hl the laser interaction between atomic or molecular 
states. From the time-independent Schrodinger equation 
e = E^J e under Born-Oppenheimer approximation, 
one obtains the following coupled equations 
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where is the laser frequency, Ai and A2 are the laser- 
induced transition dipole matrix elements between | e) 
and I gx), and between | e) and | 32), respectively. Here 
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FIG. 2: PA rate Kpa in cm 3 sec 1 is plotted as a function of 
magnetic field B in Gauss. The parameters are T p = 0.1 MHz 
(solid line), T p = 0.5 MHz (dashed line) and F p = 1.0 MHz 
(dashed-dotted line). The parameter is chosen to be - 8.90 
(solid line), -7.63 (dashed) and - 6.36 (dashed-dotted) so that 
the minimum and maximum of Kpa occur at B = 710 G and 
B ~ 736 G, respectively. The other parameters are 7 = 11.7 
MHz, S pc = - 2.8 MHz and temperature T = 10fjK. The laser 
is assumed to be tuned on resonance with the continuum- 
bound transition. 



Vi ( i = 1,2) are the potentials including hyperfine and 
Zeeman terms, V e is the excited state molecular poten- 
tial and V stands for hyperfine spin coupling between 
closed channel bound state and continuum states. Here 
Bj(r) = h 2 J(J+ l)/(2^ir 2 )] is the rotational term of the 
excited state. Note that for the ground scattering and 
bound states we have considered only the zero rotational 
state. The zero of energy scale is taken to be the thresh- 
old Eth of the open channel 1 and the energies of the 
bound states are measured from this reference. For two 
homonuclear atoms, the asymptotic form of the potential 
V e (r — » 00) ~ huoA—Cz/r^, where C3 is the long-range co- 
efficient of dipole-dipole interaction between one ground 
state S-atom and another excited state P-atom and lua 
is the atomic frequency. These three coupled equations 
can be solved exactly by the use of real space Green's 
function as described below. 

It is convenient to write V e = hujA+V ex , where V ex (r — > 
00) ~ —Cz/r 3 . Let (j) p denote the bound state solution 
of the the potential V ex and E p be the corresponding 
bound state (negative) energy. The Green's function for 
homogeneous part with Ai = A2 = (i.e. without laser 
couplings) of Eq. ([2]) can be written as 



G p {r,r') = - 
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HS + E - E p + ihj/2 



MrW) (5) 



where 5 — lol ~ loa- Using this function, we can write 
down the solution of equation @ in the form $ p = 
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A p 4>p(r) where 

fdr> [A 1 (r')x(r-')+A 2 (r')^(r')]^(/) 



-4, 



M + E - E p + ih-y/2 



(6) 



Similarly, with the use of Green's function for the ho- 
mogeneous part of Eq. ([3|), we have <&/(?") = A m <j)f(r) 
where 

, = J dr' lV*(r>) X (r>) + A p A 2 (/)0 p (r')] <^(r') 

E-Ef [ ' 

where 4>f(r) is the wave function and Ef is the energy 
of bound state in the closed channel in the absence of 
laser field. Now, we can express A p in terms of integrals 
involving the continuum state x an d molecular bound 
states 4> p and <j)f. Then substituting $ p and &f expressed 
in terms of A p and \ m Eq. ^ and making use of the 
relation \ = J dEbEipE, we obtain 

-f;J^(r) + [Vi{r)-E]i> E {r) = -ht{r)A p 4> p {r) 



{V fc + A p A p f) 
E-Ef 



V(r)0 p (r) 



where 



A pc (E - E f ) + Ap f V fc 
V(E-E f )-\A pf \* 



(8) 



(9) 



Here V = US + E — E p + ihrf/2, V fc = 
J dr(pf(r)V(r)tpE{r), A pc = J dr^ p (r)Ai(r)^ E (r) and 
A-2,pf — J dr<fip(r)A2(r)<fif(r). Equation (J8j> can now 
be solved by constructing the Green's function with the 
scattering solutions of the homogeneous part (i.e., for 
Ai = V — 0). This Green's function can be written as 

K(r,r') = -^T 9 {r)^ rr (.r>) + H> / eg (r)^/ eg (r% 

(r' > r) 

K(r,r>) = -^T 9 ir')^ E irr {r) + i^ 9 (r)^ / e9 (/)], 

(r' < r) 

where the regular function ^ E re9 (r) vanishes at r = and 
the irregular solution ip E ^ rr (r) is defined by boundary 
only at r — > oo. These have the familiar asvmptotic 
behavior ip E re9 (r) ~ jo cos 770 — no sin 770 an d ij) E rr (r) ~ 
no cos 770 + jo sin 770 , where jo and no are the spherical 
Bessel and Neumann functions for I — and 770 is the 
s-wave phase shift in the absence of laser and magnetic 
field couplings. Here E — h 2 k 2 /(2fi) with /j, being the 
reduced mass of the two atoms. Next, we can express 
the solution of Eq. ([5]) in the following form 



4> E = exp(i m yj° E re9 + I dr'K{r,r') A$(r')A p <p p (r') 
Vfc + ApApf 



E-Ef 



-V(r')Mr') 



(10) 



The stimulated line width of photoassociated 
molecule is given by the Fermi- Golden rule ex- 
pression T p — 2ir\A p ' c \ 2 /h and the Feshbach res- 



onance line width is T 



J 



2n 

0,reg , 



fc 



/h, 



Jdr(t> f (r)V(r)^ re9 (r) and A° 



where V® c 

Jdr</)p(r)Ai(r)^ re9 (r). The Stark energy shift due to 
laser coupling of PM state with the continuum is given 

by S pc = J J dr'dr^ p (r)AJ(r)Re[X;(r / ,r)]Ai(r')0p(r / )- 
Further, the physics of Feshbach resonance leads us to 
introduce the parameter 



v Pf 



dr 1 dr4>f(r)V(r)Re[K.{r' ,r)]Al(r')(t> p {r') 



which represents an effective continuum-mediated 
magneto-optical coupling between the two bound states 
where S fc = JJdr'dr^ f (r)V*(r)Re[K(r',r)]V(r')^ f (r') 
is the energy shift of the closed channel bound state 
due to its coupling with the continuum. Now writing 
e = (E - Ef)/(Y } /2) with E f = E f + S fc being the 
shifted energy of the closed-channel bound state, and in- 
troducing a parameter 



A 



pf 



Vr 



pf 



,u T)c v fc 



(11) 



which we call "Feshbach asymmetry parameter" , we can 
express 



, = ( e + q f ^ cxp(777 ) ^ 



HTf/2 



Ap+i-f/Tf+D! 



[E - (E p - hS)]/(T f /2) 



[E P 



where A p 

HS — Ef]/(T f /2) is independent of laser intensity and 
Dj = {-2S pc + %T p )/r f - (T p /T f ){q f - i) 2 /{e + i) is a 
parameter which is proportional to laser intensity /. In 
writing the above equation we have assumed that V^ c , 

and Ap C are real quantities. Note that q/ is independent 
of laser power since its numerator as well as the denom- 
inator is proportional to laser amplitude. Following the 
Ref. [23| , we can express e in terms of applied magnetic 
field in the form 

e _ E - E t h — (E t h - Ef) _ E - E th _ B — B 



17/2 



Tf/2 A(ka bg ) 



where E t h is the threshold of the open channel, A is the 
Feshbach resonance width, Bq is the resonance magnetic 
field and a bg is the background scattering length. Here 
E — E t h is the asymptotic collision energy. The energy 
E t h depends on the applied magnetic field due to Zeeman 
shift of the atomic level. The resonance scattering length 
is given by a res = -a bg A/ (B - B ). 

Before we discuss our main results, we would like to 
point out how our mathematical treatment discussed 
above is related to the recent work of Koznetsova et 
al. [2(| who have studied a related model in a differ- 
ent context which is to transfer of atoms into ground 
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FIG. 3: The dashed lines are Kpa vs. B plots when PA laser 
is tuned on resonance with continuum-bound transition. The 
solid lines are corresponding plots when the laser is tuned 
on resonance with bound-bound transition. For upper panel 
(a) T p = 0.01F/, and for lower panel (b) F p = 0.05F/. For 
both the panels, the other fixed parameters are g/ = —9.0, 
7 = 0.7T f , T = lOuK, T f = 16.67 MHz and B = 736.72 
Gauss. The minimum and maximum of dashed curve in (a) 
appear at 710 Gauss and 736.8 Gauss, respectively. 



state molecules via two-photon process near a Feshbach 
resonance. Our approach is to find out the real space 
dressed wave function by solving time-independent scat- 
tering problem by Green's function method while they 
have adapted a quantum optics-based approach of find- 
ing time-dependent amplitudes of the dressed state by 
solving coupled differential equations numerically. 



3. THE RESULTS AND DISCUSSION 

We now discuss characteristic features of our main re- 
sults. PA spectrum is given by the PA rate coefficient 
Kp A = {vreiopA}, where a PA = (hj\A p \ 2 )/(2TTk 2 ) is the 
cross section for the loss of atoms due to decay of the 
excited molecules. Here (• • •} implies thermal averaging 
over the relative velocity v re \ = hk/fi. Note that, in the 
limit Tf — ► 0, PA spectrum reduces to a Lorentzian im- 
plying that coupling between closed channel bound state 
and the continuum is essential for the occurrence of Fano 
interference. When both A2 and Vf c go to zero, the spec- 
trum reduces to that of standard PA. For numerical illus- 
trations, we consider a model system of two ground-state 
(S1/2) 7 Li atoms undergoing PA from the ground molec- 
ular configuration 3 E+ to the vibrational state v = 83 
of the excited molecular configuration 3 S g which corre- 
lates asymptotically to 2S1/2 + ^1/2 f ree atoms. The 
spontaneous linewidth is taken to be 7 = 11.7 MHz 
24]. The experimental value of shift S pc is reported be 
-1.7 ± 0.2 MHz /W cm 2 The resonance width is 

A = —192.3 Gauss and the background scattering length 
dbg = — 24.5ao (ao is Bohr radius) [2||. The Feshbach 



resonance linewidth Tf at 10 /iK temperature is calcu- 
lated out to be 16.66 MHz using the parameters reported 
in Ref. 

Depending on how PA laser is tuned, we have two 
cases. In the first case (case-I), laser is on or near res- 
onance with free-bound transition but off-resonant with 
bound-bound transition. In the second case (case-II), it 
is resonant with bound-bound transition. PA rate will 
be maximized at the poles of Eq. (fT!?]) . In case-I, there is 
only one pole of Eq. (I12p which depends on laser inten- 
sity. The minimum in the spectrum is solely determined 
by the asymmetry parameter g/ and is independent of 
laser intensity. We first consider case-I and plot (Kpa) 
as a function of B for three different values of T p in Fig. 2. 
For these three different T p values we choose three differ- 
ent qf parameters such that the maximum appears near 
B = 736 G [l4j . Since the minimum position B m i n is in- 
dependent of laser intensity, we also choose three different 
values of resonant magnetic field Bq such that B m i n re- 
mains fixed at 710 G for theses three T p values. In this 
case there arises asymmetric Fano profile with one mini- 
mum and one maximum. This results from quantum in- 
terference between continuum-bound and bound-bound 
Raman-type transition pathways. This interpretation of 
Raman Fano profile is in accordance with the recent ex- 
perimental observation of two-photon PA by Moal et al. 

Next we consider case-II in which PA laser is tuned 
in resonance with bound-bound rather than continuum- 
bound transition. In this case we have HS — E p + Ef = 
and so A p = e. Then A p will have two maxima given by 

e(e + i) + if *(e + i) - f p (q f - if = 0. (14) 

where T t = T t /Tf with T t = T p + 7 being the total 
line width, and T p = T p /Tf. For the sake of compari- 
son, we plot spectra in Fig. 3 for both the cases. Figure 
3(a) shows that a single maximum appears in both the 
cases when laser intensities are low. As laser intensity 
increases, the maximum in case-I disappears while a two 
peak structure emerges in case-II as displayed in Fig. 3(b). 
We notice that PA rate is lower in case-II in comparison 
to case-I for the same magnetic field and other parame- 
ters except near the two maxima. To further investigate 
into the double-peak structure, we demonstrate spectra 
for case-II at higher intensities in Fig. 4 which clearly in- 
dicates the nonlinear features of Fano interference. The 
origin of the two peaks lies in Autler-Townes splitting Q 
due to bound-bound resonant coupling when continuum- 
bound dipole coupling is scanned into two-photon reso- 
nance. At lower intensities, the two peaks can appear 
on the same side of Fano minimum. As a result, there 
can appear another smaller minimum (we call it Autler- 
Townes (AU) minimum in order to distinguish it from 
Fano minimum) between the two peaks. By comparing 
Fig. 3 (a) with Fig. 3(b), we note that the AU minimum at 
a higher intensity appears near the position where max- 
imum would have appeared at a lower intensity. The 
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FIG. 4: Shown are Kpa vs. B plots when PA laser is tuned 
on resonance with bound-bound transition for F p = 0.57 ( a )i 
r p = 7 (b), Tp = 27 (c) and F p — 47 (d). The parameters 
chosen are Bq = 730.51Gauss, qf — —6.89 with all others 
parameters remaining same as in Fig. 2 

separation between the two maxima increases with in- 
creasing laser intensities as shown in Fig. 4. As one of the 
peaks crosses the Fano minimum at an increased laser 
intensity, the AU minimum disappears due to its inter- 
ference with the much stronger Fano minimum resulting 
in two-maximum structure only. The double-maximum 
structure is particularly prominent in the strong-coupling 
regime where T p exceeds the spontaneous line width 7 
of PA molecule. Recently, Pellegrini and Cote [18| have 
theoretically obtained double-minimum spectra using the 
formalism of Ref. which is to first diagonalize the 
part of the Hamiltonian pertaining to the ground state 
scattering (continuum interacting with the bound state 
in closed-channel), and then to calculate the optical tran- 



sition matrix element between this diagonalized state and 
the excited molecular state by Fermi- Golden rule. This 
is similar to linear Fano theory [l[ and hence can not be 
applied for strong-coupling that can further modify the 
continuum state significantly. In our formalism, we have 
diagonalized the full Hamiltonian nonpcrturbatively. 



4. CONCLUSION 

The results discussed above clearly demonstrate linear 
and nonlinear aspects of Fano interference in weak- and 
strong-coupling regimes, respectively. Observation of 
two-minimum and two-maximum structures crucially de- 
pends on precise tuning of PA laser on or near resonance 
with bound-bound transition. If the laser field is tuned 
to get the maximum amount of loss of atoms for a fixed 
magnetic field, the resulting PA spectrum will mostly 
correspond to the case-I with a single maximum. To 
explore the nonlinear Fano effect, it is important to know 
the binding energy of the closed-channel bound state so 
that the laser can be accurately tuned near resonance 
with the bound-bound transition as the magnetic field is 
varied. Recently, nonlinear Fano effect was observed in 
quantum dot [27| . Although Autler-Townes split ting has 
been recently demonstrated in two-photon PA [1JJ, [UJ , 
it is yet to be observed in PA with a single laser beam 
in the presence of Feshbach resonance. Fano interference 
may further be explored in photoassociation between 
heteronuclear atoms such as Na and Cs [28] or K and 
Rb [HI which have broad magnetic Feshbach resonance 
and shorter ranged excited potentials. 
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